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1.0  INTRODUCTION 

The  computation  of  reentry  vehicle  motion  requires  knowledge  of  a 
number  of  aerodynamic  coefficients;  these  may  come  from  correlations  of 
experimental  data  or  from  theoretical  calculations.  As  more  accurate 
trajectory  predictions  are  required,  more  accurate  aerodynamic  coeffi- 
cients must  be  obtained.  Because  many  vehicles  are  caused  to  spin 
during  entry  into  the  atmosphere,  the  effects  of  viscous  aerodynamic 
damping  on  the  spin  rate  must  be  known  in  order  to  predict  the  per- 
formance of  the  vehicles,  A knowledge  of  roll-damping  effects  is  also 
valuable  for  the  prediction,  corroboration,  and  evaluation  of  wind 
tunnel  test  data. 

Approximate  methods,  such  as  those  used  in  Ref.  1,  have  been  applied 
in  the  past  to  compute  roll-damping  data,  and  such  methods  do  provide 
reasonable  results  when  properly  applied.  Recently,  more  exact  methods, 
such  as  in  Ref.  2 and  those  used  in  the  present  study,  have  become 
available  for  computation  of  the  viscous  flow  over  spinning  bodies. 

It  is  the  purpose  of  this  report  to  present  a method  which  has  been 
developed  for  the  computation  of  the  boundary- layer  flow  over  spinning 
sharp  and  blunt  cones  in  supersonic  and  hypersonic  streams,  and  to 
present  the  results  of  a parametric  study  of  the  effects  of  Mach  number, 
Reynolds  number,  cone  angle,  and  bluntness  ratio  on  the  roll-damping 
derivative  for  a family  of  sharp  and  blunt  cones  at  zero  angle  of  attack. 
Both  laminar  and  turbulent  boundary  layers  have  been  considered,  and  for 
the  blunt  cases  the  effects  of  swallowing  of  the  inviscid  entropy  layer 
by  the  boundary  layer  were  included  in  the  analysis.  Some  investigations 
were  also  made  of  the  effects  of  vehicle  wall  temperature  and  roll  rate, 
and  Reynolds  number  scaling  of  the  roll-damping  derivative  was  examined. 
Conditions  were  considered  which  are  typical  of  both  wind  tunnel  and 
free-f light  situations,  and  limited  comparisons  were  made  with  experi- 
mental data  and  results  of  other  computations. 
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2.0  THEORETICAL  METHOD 

The  treatment  of  the  boundary- layer  equations  used  in  this  investi- 
gation evolved  from  a study  of  three-dimensional  boundary- layer  flows 
(Refs.  3 and  4).  This  resulted  from  the  fact  that  the  three-dimensional 
boundary-layer  equations  include  the  lateral  momentum  equation  and 
contain  the  cross-flow  velocity  terms  needed  to  treat  the  spinning  body 
case,  although  they  are  much  more  general  than  needed  in  the  axisymmetric 
flow  situations  considered  in  this  investigation. 

The  compressible  boundary- layer  equations,  including  the  cross-flow 
velocity  terms,  were  cast  into  a normalized  Crocco-variables  form  and 
solved  using  an  implicit  finite  difference  technique  developed  by  Taylor 
series  truncation.  The  Crocco  transformation  is  accomplished  by  letting 
the  coordinate  normal  to  the  body  surface  be  the  normalized  longitudinal 
velocity  component.  This  yields  the  convenient  limits  of  zero  and  unity 
for  the  normal  coordinate.  The  velocity  component  normal  to  the  body 
surface  is  eliminated  from  the  governing  equations  by  combining  the 
longitudinal  momentum  equation  with  each  of  the  continuity,  lateral 
momentum,  and  energy  equations.  This  decreases  the  number  of  partial 
differential  equations  to  be  solved  from  four  to  three.  As  a part  of 
the  transformation,  a form  of  the  longitudinal  component  of  the  shearing 
stress  becomes  the  dependent  variable  associated  with  the  transformed 
longitudinal  momentum/continuity  equation.  The  energy  equation  is  cast 
in  terms  of  the  total  enthalpy,  and  the  gas  under  consideration  has  been 
treated  as  both  thermally  and  calorically  perfect,  being  air  for  all  of 
the  cases  considered.  The  specific  forms  of  the  governing  equations, 
transport  laws,  difference  equations,  etc.,  are  not  given  herein,  but 
are  the  axisymmetric  flow  forms  of  those  given  in  Refs.  3 and  4. 

For  turbulent  boundary-layer  flow,  the  viscosity  and  Prandtl  number 
in  the  governing  equations  were  treated  as  effective  transport  parameters, 
and  an  eddy-exchange  coefficient  treatment  of  the  turbulent  transport 
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parameters  was  used.  The  turbulent  viscosity  was  specified  using  the 
mixing-length,  invariant  turbulence  model  as  applied  by  Hunt,  Bushnell, 
and  Beckwith  (Ref.  5),  with  the  mixing  length  obtained  using  a two- layer 
model  with  exponential  damping  near  the  wall  as  recommended  by  van 
Driest  (Ref.  6).  The  variation  of  the  mixing  length  was  obtained  from 
the  recommendations  of  Patankar  and  Spalding  (Ref.  7).  Transition  from 
laminar  to  turbulent  flow  was  treated  by  the  use  of  the  streamwise 
intermittency  factor  determined  by  Dhawan  and  Narashima  (Ref.  8),  with 
the  beginning  and  end  of  transition  being  fixed  at  x/r^  = 1.5  and  2.0, 
respectively,  for  the  blunt  cases  and  at  x = 0. 1 ft  and  0.15  ft,  respec- 
tively, for  the  sharp  cases.  A turbulent  Prandtl  number  of  0.9  and  a 
laminar  Prandtl  number  of  0. 7 were  assumed  in  this  work. 

The  governing  equations  degenerated  to  two-point  boundary-value 
problems  at  a sharp  tip  or  stagnation  point,  and  formed  a parabolic 
system  which  could  be  .solved  by  marching  down  the  remainder  of  the  body. 

The  boundary  conditions  required  were  generally  the  values  of  the  depen- 
dent variables  on  the  body  surface  and  at  the  outer  edge  of  the  boundary- 
layer,  together  with  the  longitudinal  derivatives  of  the  outer-edge 
conditions.  Wall  boundary  conditions  were  provided  by  the  no-slip 
condition  and  a specified  surface  temperature,  together  with  the  in- 
viscid  surface  pressure,  which  was  assumed  constant  with  respect  to  the 
direction  normal  to  the  surface,  the  standard  boundary- layer  assumption. 

For  the  sharp  cases,  the  outer-edge  boundary  conditions  were  taken 
as  the  inviscid  body  surface  conditions.  For  the  blunt  cases,  the 
outer-edge  boundary  conditions  were  determined  from  inviscid  flow  data, 
based  on  consideration  of  the  swallowing  by  the  boundary  layer  of  the 
inviscid  entropy  layer.  Entropy  swallowing  is  based  on  a mass-flow  balance 
between  the  boundary  layer  at  a given  location  and  the  free  stream  in 
order  to  find  the  shock  angle  crossed  by  the  outer-edge  flow  at  the 
location  of  interest.  This  is  shown  in  Fig.  1.  Inviscid  shock  shape 
and  surface  pressure  data  were  obtained  using  the  blunt  body  treatment 
of  Ref.  9 and  the  method  of  characteristics  treatment  of  Ref,  10. 
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Figure  1.  Inviscid  entropy  layer  swallowing  by  boundary  layer. 


This  method  of  solution  of  the  governing  partial-different  equa- 
tions involved  replacing  them  by  a set  of  consistent,  linearized,  alge- 
braic equations.  The  resulting  set  was  of  tridiagonal  form  and  could  be 
solved  by  means  of  the  well-known  algorithm  available  for  such  cases. 

The  solution  required  iteration  to  remove  the  assumptions  made  in  the 
linearization,  and  cases  which  considered  entropy  swallowing  involved  an 
outer  iteration  concerned  with  the  determination  of  consistent  outer- 
edge  boundary  conditions. 


3.0  RESULTS  OF  CALCULATIONS 

In  this  section,  results  of  calculations  made  using  the  method 
described  in  the  previous  section  are  presented.  Comparisons  are  made 
with  experimental  data  and  results  from  another  method  of  calculation. 
These  comparisons  tend  to  validate  the  method  developed  in  this  inves- 
tigation. Results  of  an  examination  of  the  effects  of  vehicle  wall 
temperature  and  spin  rate  on  the  roll-damping  derivative  are  presented, 
and  some  considerations  of  Reynolds  number  scaling  are  also  discussed. 
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The  bulk  of  the  results  obtained  in  this  investigation  are  pre- 
sented as  a series  of  figures  which  give  the  behavior  of  the  roll- 
damping  derivative  of  a sphere-cone  as  a function  of  body  bluntness 
(ratio  of  nose  radius  to  base  radius)  for  a range  of  Mach  -numbers, 

Reynolds  number,  and  cone  angles.  Results  are  shown  for  both  laminar 
and  turbulent  boundary- layer  flows  for  those  cases  pertinent  to  wind 
tunnel  flows,  namely  Mach  numbers  of  6,  8,  10,  and  14.  For  the  Mach  22 
data  which  is  representative  of  atmospheric  entry  conditions,  laminar 
flow  results  are  given  for  altitudes  of  200,  120,  and  80  kft,  and  turbulent 
flow  results  are  given  for  altitudes  of  120,  80,  and  50  kft. 

All  of  the  calculations  made  in  this  investigation  were  performed 
on  an  IBM  370/165  digital  computer. 

3.1  COMPARISONS  WITH  OTHER  DATA 

To  validate  the  method  used,  results  of  calculations  made  using  the 
present  method  have  been  compared  with  experimental  data  and  with  results 
from  another  method  of  computation.  Figure  2 shows  the  roll-damping 
derivative  versus  nose  bluntness  ratio  for  laminar  flow  over  a rotating 
10-deg  blunt  cone  in  a Mach  14.25  flow.  Results  from  boundary-layer 
calculations  are  compared  with  experimental  data  reported  in  Ref.  1. 

The  agreement  between  the  measurements  and  the  calculations  is  reason- 
ably good,  with  the  computed  results  being  from  one  to  five  percent 
above  the  measurements,  depending  on  the  bluntness  ratio.  It  is  dif- 
ficult to  make  a statement  about  the  quantitative  agreement  between  the 

* 

two  sets  of  data  because  of  a lack  of  information  about  the  uncertainty 
in  the  experimental  data.  It  is  known  that  accurate  experimental  data 
of  this  nature  are  hard  to  obtain  because  of  the  large  beating-damping 
corrections  which  must  usually  be  made.  The  present  method  was  further 
validated  by  solving  the  sharp-cone  case  of  Fig.  2 using  the  method  of 
Ref.  2.  The  results  obtained  using  that  method  agreed  to  within  two 
percent  with  the  lateral  skin-friction  coefficient  obtained  using  the 
present  method. 
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Figure  2.  Measured  and  computed  roll-damping  derivative  versus 

bluntness  ratio  for  a 10-deg  cone  in  a Mach  14.25  flow. 

In  general,  the  calculations  made  in  this  investigation  used  the 
Sutherland  viscosity- temperature  relationship;  however,  for  the  cases 
shown  in  Fig.  2,  the  Sutherland  law  was  used  only  for  temperatures  above 
200°R,  with  a linear  temperature-viscosity  relationship  being  used  for 
temperatures  below  200°R.  This  was  necessitated  by  the  fact  that  Che 
free-stream  static  temperature  for  these  cases  was  approximately  50°R, 
and  the  viscosity-temperature  relationship  is  not  well  modeled  at  this 
low  temperature  by  the  Sutherland  Law.  Since  the  free-stream  conditions 
are  used  as  reference  conditions  in  the  theoretical  treatment  of  the 
problem  it  is  necessary  that  the  free-stream  viscosity  be  properly 
computed  even  though  such  low  temperatures  do  not  actually  exist  within 
the  boundary  layer.  Use  of  the  Sutherland  viscosity  law  exclusively  for 
the  cases  of  Fig.  2 yielded  a roll-damping  derivative  curve  approx- 
imately thirteen  percent  above  that  shown.  For  the  remainder  of  the 

i 

calculations  made  in  this  investigation,  the  free-stream  static  temper- 
ature was  100°R  or  higher,  and  for  these  cases  only  the  Sutherland 
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viscosity  was  used.  The  estimated  variation  in  the  roll- damping  de- 
rivative is  less  than  three  precent  when  ignoring  the  linear  portion  of 
the  viscosity-temperature  law  for  these  cases. 

One  further  point  to  be  made  in  examining  Fig.  2 is  the  manner  in 
which  the  blunt  cone  results  approach  the  sharp  cone  (r^/r^  = 0)  value 
as  the  bluntness  ratio  is  reduced.  Physically,  this  is  a natural  occur- 
rence, but  computationally  it  occurs  only  because  of  the  inclusion  of 
entropy  swallowing  in  the  analysis.  In  fact,  if  entropy  swallowing  were 
neglected,  the  blunt-cone  results  would  approach  a value  approximately 
fifteen  percent  below  the  sharp-cone  value  as  the  bluntness  is  reduced 
to  zero.  Since  the  sharp-cone  results  were  obtained  using  an  entirely 
different  code  than  for  the  blunt  cone,  the  consistency  between  the 
sharp-cone  value  and  the  results  obtained  including  entropy  swallowing 
in  the  blunt-body  analysis  demonstrates  a consistency  between  the  two 
codes  and  the  necessity  for  including  entropy  swallowing  in  the  compu- 
ations. 

3.2  BODY  GEOMETRY  AND  FREE-STREAM  CONDITIONS 

The  remainder  of  the  results  to  be  presented  are  roll-damping 
derivative  data  on  a family  of  sharp  and  blunt  cones  of  semi-vertex 
angles  of  5.0,  7.5,  and  10.0  deg.  For  a given  cone  angle,  the  base 
radius  was  held  constant,  and  the  nose  radius  varied  to  give  various 
bluntness  ratios.  Base  radii  of  4.0,  4.5,  and  5.0  in.  were  used  for  the 
5.0-,  7,5-,  and  10.0-deg  cones,  respectively.  Calculations  were  per- 
formed for  nose-to-base  radius  ratios  of  0.0,  0.05,  0.15,  and  0.25  in 
each  case. 

For  the  Mach  6,  8,  10,  and  14  cases,  a free-stream  static  tem- 
perature of  100°R  and  unit  Reynolds  numbers  of  0.5  x 10^  and  3.0  x 10^ 
per  ft  were  considered  as  being  typical  of  wind  tunnel  conditions. 
Standard  values  of  wall  temperature  of  600°R  and  spin  rate  of  500  rpm 
were  used.  Laminar  and  turbulent  flow  were  considered  in  all  of  these 
cases. 
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For  the  free-flight  cases,  a Mach  number  of  22  was  used  together 
with  a wall  temperature  of  2000°R  and  a spin  rate  of  100  rpm.  The  same 
geometries  as  previously  mentioned  were  treated.  Cases  were  run  for 
laminar  flow  at  200-,  120-,  and  80-kft  altitudes,  and  for  turbulent  flow 
at  altitudes  of  120,  80,  and  50  kft. 

3.3  WALL  TEMPERATURE  AND  SPIN-RATE  EFFECTS 

Variations  of  wall  temperature  and  spin  rate  from  that  of  Section 
3.2  were  considered,  and  the  results  are  presented  in  this  section. 

The  Mach  8,  = 7.5-deg,  r^/r^  =0.15  case  was  run  with  spin  rates 

of  100,  500,  and  1,000  rpm,  and  with  wall  temperatures  of  600°R  and  an 
adiabatic  wall  condition.  These  cases  were  run  at  unit  Reynolds  numbers 
of  0.5  x 10  and  3.0  x 10  /ft,  considering  both  laminar  and  turbulent 
boundary  layers. 

For  these  various  conditions,  the  net  effect  of  the  spin  rate 
variation  on  the  roll-damping  derivative  was  less  than  one  percent  for 
the  range  of  spin  rate  considered,  holding  other  parameters  fixed. 

The  roll-damping  derivative  was  substantially  increased  by  changing 
the  wall  temperature  from  600°R  to  an  adiabatic  wall  condition  (a  wall 
temperature  in  the  vicinity  of  1200°R).  For  the  lower  Reynolds  number 
laminar  and  turbulent  cases  and  the  higher  Reynolds  number  laminar 
cases,  this  increase  was  approximately  ten  percent.  For  the  higher 
Reynolds  number  turbulent  cases,  the  increase  caused  by  the  increase  in 
wall  temperature  was  approximately  twenty  percent. 

3.4  REYNOLDS  NUMBER  SCALING 

Typical  results  of  the  parametric  study  of  the  effects  of  Mach 
number,  Reynolds  number,  cone  angle,  and  bluntness  ratio  on  the  roll- 
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damping  derivative  of  blunt  cones  are  shown  in  Fig.  3.  This  figure 
shows  the  roll-damping  derivative  versus  nose  bluntness  ratio  for  both 
laminar  and  turbulent  flow  on  a 7.5-deg  cone  in  a Mach  8 free  stream. 
Because  of  (a)  the  inclusion  of  the  cross-flow  terms  in  the  governing 
equations,  (b)  the  consideration  of  entropy  layer  swallowing  by  the 
boundary  layer,  and  (c)  the  consideration  of  turbulent  flow,  the  governing 
equations  considered  in  this  study  do  not,  strictly  speaking,  have 
similarity  solutions  which  allow  casting  the  results  in  a Reynolds 
number  independent  form.  Evidently,  however,  the  effects  of  (a)  and  (b) 
are  not  strong  in  the  laminar  cases  shown  in  Fig.  3,  and  the  application 
of  conventional  Reynolds  number  scaling  yields  the  essentially  Reynolds 
number  independent  data  shown  in  the  upper  part  of  Fig.  4.  Heuristically 


Figure  3.  Rol 1 -damping  derivative  versus  bluntness  ratio  for  a 
7.5-deg  blunt  cone  in  a Hach  8 stream  under  laminar 
and  turbulent  flow  conditions. 

treating  the  turbulent  flow  data  of  Fig.  3 in  a similar  manner  yields 
the  result  shown  in  the  lower  part  of  Fig.  4,  where  the  0.27  power  on 
the  Reynolds  number  was  determined  empirically  as  a value  which  reason- 
ably collapsed  the  data.  Although  the  remainder  of  the  data  to  be 
presented  in  this  report  will  be  in  the  format  of  Fig.  3,  having  the 
data  in  the  form  given  in  Fig.  4 could  make  its  use  more  convenient  in 
motion  analysis. 
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Figure  4.  Roll-damping  derivative,  scaled  by  Reynolds  number  to  a 
power,  versus  bluntness  ratio. 

3.5  PARAMETRIC  STUDY  RESULTS 

The  bulk  of  the  results  obtained  in  this  investigation  are  pre- 
sented in  Figs.  5 through  10.  These  figures  give  the  variation  of  the 
roll-damping  derivative  of  a sphere-cone  as  a function  of  body  blunt- 
ness. Each  figure  has  a separate  set  of  curves  for  cone  angles  of  5.0, 

7.5,  and  10.0  deg. 

Figures  5 through  8 present  data  for  f ree-s.tream  Mach  numbers  of  6, 

8,  10,  and  14,  respectively.  Each  Mach  number-cone  angle  combination 

has  both  laminar  and  turbulent  flow  results  for  free-stream  unit  Reynolds 

6 6 

numbers  of  0.5  x 10  and  3.0  x 10  /ft.  Generally  speaking,  the  negative  of 
the  roll-damping  derivative  decreases  with  increasing  bluntness.  The  varia- 
tion between  the  sharp  and  the  0.25  bluntness  values  is  generally  greater 
for  the  turbulent  cases  than  the  laminar,  i.e.,  the  effect  of  bluntness  is 
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greater  for  the  turbulent  cases.  Also,  for  the  range  of  bluntness 
considered,  the  effect  of  bluntness  generally  increases  with  increasing 
free-stream  Mach  number. 


rn'rb 


c.  0 = 10.0  deg 

w 

Figure  5.  Roll-damping  derivative  versus  bluntness  for 
5.0-,  7.5-,  and  10.0-deg  cones  in  Mach  6 flow. 
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Figures  9 and  10  present  the.  roll-damping  derivative  results  obtained 
for  the  Mach  22  flight  conditions  considered.  For  these  cases,  the 
free-stream  Reynolds  number  and  static  temperature  depend  upon  the 
altitude  for  each  specific  case.  The  free-stream  unit  Reynolds  number 
varies  from  approximately  26  million  per  ft  at  50-kft  altitude  to  36 
thousand  per  ft  at  200-kft  altitude,  and  the  mean  free-stream  tempera- 
ture is  approximately  420°R,  Figure  9 presents  the  results  for  the  50-, 
80-,  and  120-kft  altitude  cases,  where  both  laminar  and  turbulent  results 
are  given  for  the  two  higher  altitudes,  and  only  turbulent  results  are 
given  for  the  lowest  altitude  case.  In  Fig.  10,  laminar  flow  results 
are  presented  for  the  200-kft  altitude  cases. 

The  primary  purpose  of  Figs,  5 through  10  is  to  provide  data  for 
use  in  the  analysis  of  the  motion  of  spinning  reentry  vehicles  in  the 
wind  tunnel  and  in  free  flight.  It  is  hoped  that  the  range  of  param- 
eters considered  is  adequate  for  this  purpose. 

4.0  SUMMARY  AND  CONCLUSIONS 

A theoretical  method  of  computing  the  boundary- layer  flow  over 
spinning  sharp  and  blunt  cones  in  supersonic  and  hypersonic  flow  has 
been  described.  The  method  was  validated  by  being  shown  to  give  good 
agreement  with  experimental  data  and  another  method  of  computation. 

Wall  temperature  and  spin-rate  effects  on  roll-damping  derivative  data 
were  examined,  and  a method  of  casting  the  roll-damping  derivative  data 
into  a Reynolds  number  independent  form  was  presented. 

The  bulk  of  the  data  presented  were  the  results  of  a parametric 
study  of  the  effects  on  the  roll-damping  derivative  of  Mach  number, 
Reynolds  number,  cone  angle,  and  bluntness  ratio.  The  results  were  for 
a range  of  the  parameters  which  should  allow  their  application  to  the 
analysis  of  the  motion  of  spinning  cones  in  both  wind  tunnel  and  flight 
situations. 
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Figure  10.  Ro 11 - damp i ng  derivative  versus  bluntness  for  5.0-, 
7.5-,  and  10.0-deg  cones  in  Mach  22  flow  at  200-kft 
altitude  with  laminar  flow. 
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NOMENCLATURE 


Base  area 

Roll-damping  derivative, 


Base  diameter 

Rolling  moment 

Free-stream  Mach  number 

Roll  rate,  radians/sec 

Free-stream  dynamic  pressure 

Free-stream  unit  Reynolds  number,  1/ft 

Free-stream  Reynolds  number  based  on  base  radius 

Base  radius 

Nose  radius 

Ratio  of  surface  temperature  to  free-stream  total 
temperature 

Free-stream  velocity 

Surface  distance  from  stagnation  point  or  sharp 
tip 

Cone  angle  (semi-vertex  angle) 

Roll  rate,  rev/min 
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